Raman resonance in spin S two-leg ladder systems 
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We argue that the Raman intensity in a spin 5* two-leg spin-ladder has a pseudo-resonance peak, 
whose width is very small at large S. The pseudo-resonance originates from the existence of a local 
minimum in the magnon excitation spectrum, and is located slightly below twice the magnon energy 
at the minimum. The physics behind the peak is similar to the excitonic scenario for the neutron 
and Raman resonances in a d— wave superconductor. 
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Recently, there has been a considerable experimental 
progress in Raman studies of two-leg spin- ladder materi- 
als {Sr,La) li Ca x Cu 2 iO il (Ref 0,12) and SrCu 2 3 0. 
The most intriguing experimental result is the discovery 
of a sharp peak in the Raman intensity of Stx^Cu-^O^x 
at a frequency near 3000cm -1 (about AOOmeV). The 
peak exists for polarizations of incoming and outgoing 
light both along and across the ladders (xx and yy, re- 
spectively), and its width for xx polarization is around 
130cm -1 which is nearly 10 times smaller than the width 
of the two-magnon peak in 2D antiferromagnets . 

The discovery of the peak stimulated the search for 
possible resonance-like features in the two-magnon Ra- 
man profile R(oj) of spin-ladder systems M El 0- 
Previous studies of S = 1/2 ladders didn't find the res- 
onance and suggested more complex explanation of the 
sharp peak in R(u) 0,0]- We argue in this paper that the 
Raman intensity in a spin S two-leg ladder possesses a 
pseudo-resonance, whose origin is similar to the origin of 
the neutron resonance and Big Raman pseudo-resonance 
in high T c superconductors 0,0- The intrinsic width of 
the Raman resonance is small in 1/5 for S > 1, but 
increases as S decreases. We discuss whether the sharp 
Raman peak observed in S = 1/2 material SV14CU24O41 
may be this resonance. 

The pseudo-resonance in a two-leg ladder emerges be- 
cause of the existence of a local minimum in the magnon 
spectrum. The magnetic excitations in a ladder are well 
described by Heisenberg interactions between spins along 
the legs of the ladder (Ji ) and along the rungs of the lad- 
der (J2). In the quasiclassical (large S) approximation, 
the excitation spectrum consists of two branches e k and 
c k +iT where 



e k = 4/1 S | sin(fc/2)| v/cos 2 (fc/2) + J 2 /2 J x 



(1) 



The spectrum e k is gapless at k = 0, and for J 2 < 2Ji, 
which we only consider, it reaches a maximum e ma x — 
S(2Ji + J2) at some finite k, then falls down and reaches 
a minimum e m in — J1J2 at k = ir (see Fig^t). At 
a finite S, Haldane effect |l2j produces a gap in e k at 
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FIG. 1: The magnon excitation spectrum in a two-leg ladder. 
Left panel, the quasiclassical case S> 1, J2/J1 = 0.4. Right 
panel - schematic e k for S — 1/2 4]. In both cases, the 
spectrum has a local minimum at k = n and a local maximum 
at some smaller k. In the quasiclassical case, the excitation 
spectrum is gapless at k = 0. For S — 1/2 ladder, magnon 
states near k — are gapped due to the Haldane effect. Note 
that our momenta are shifted by tt compared to Q- The local 
minimum disappears at J2/J1 = 2 both in the quasiclassical 
case and for S = 1/2 [jj 



k = 4], however, the minimum at k — it survives (see 
FigJIJ)). We verified that the effects on the Raman profile 
from e k and efe+ w are additive, so below we only discuss 
the dispersion branch e k We first present the summary of 
the results and then discuss the computations. Follow- 
ing previous studies of 2D systems |0, |n| , we assume 
that the Raman intensity R(ui) in two-leg ladders is rea- 
sonably well approximated by the RPA expression 0] 



R(u) = - 



ImRn 



(1 + aReRoY + (aImR y 



(2) 



where Ro(u) is the polarization bubble of free fermions 
with Raman vertices, and positive a accounts for 
magnon- magnon interaction. Near e m .; n and e m ax, exci- 
tation spectrum is flat, and the magnon density of states 
diverges. This leads to square-root singularities in Rq(uj) 
near 2c m i n and 2^. max . Near these two singularities, 

Ro(u) = 7= A i?oH ' 



V2^ 



uj - 



iS 
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FIG. 2: The behavior of the Raman bubble Ro(u>) (in units 
of Ji for non-interacting magnons. We used J2/J1 = 0.8 for 
illustrative the 
polarization bubble diverges as a square-root due to singu- 
larities in the magnon density of states. Observe that while 
ImRo(uj) is negative for all frequencies, ReRo(u>) is positive 
above 2e max , but negative below 2e min . This last behavior 
leads to a resonance in the Raman intensity below 2e m i n - In 
between 1.8e maa; and 2e max , ReRo(u>) is nonzero, but very 
small. 



FIG. 3: The theoretical behavior of the full Raman inten- 
sity R(ui) given by Eq. Left panel, the quasiclassical 
case formally extended to S = 1/2. Solid line is the result 
of numerical calculations valid at large uj, near twice mini- 
mum and maximum of the magnon spectrum. Dashed line 
is the expected behaviour for small to. We used J2/J1 = 0.8, 
as in Fig|5| Right panel, the expected behavior of R(ui) for 
S = 1/2. The width of the pseudo-resonance below 2t m i n , 
obtained in numerical simulations for S = 1/2 ladder 0,0, is 
larger than in the quasiclassical analysis extended to S — 1/2. 



where A and A* are positive. The imaginary part of 
i?o(w) diverges at approaching 2e m i n from above and at 
approaching 2e max from below (see Fig. [2J. One can 
easily make sure that in both cases ItuRq is negative 
(and Raman intensity is positive). Meanwhile, ReRo(uj) 
is positive above 2e max , and negative below 2e m i n (see 
Fig. |2J. A positive ReRo(ui) at w > 2e max implies that 
above 2e max , 1 + aReRo(uj) in Eq. is far from zero, 
i.e., there is no antibound state in the Raman profile. 
This is consistent with previous studies of the Raman 
profile in 2D Heisenberg systems ^lj. At the same time, 
below 2e max , ReRo(tu) < 0, and therefore there exists a 
frequency co res at which 1 + aReRo(cu res ) = 0, and the 
Raman intensity R(lo) is strongly enhanced. If e m i n was 
a true minimum of the excitation spectrum, the full R(u>) 
given by Eq. J2J would develop a truly S— functional res- 
onance peak at In reality, there are magnon 
states below e m i n (see FigEJ) , and ImRo remains finite, 
albeit small below 2e m i„. In this situation, the full Ra- 
man intensity acquires only a peak at lo — loq. The width 
of the peak scales as 1/S at large S, but is 0(1) for 
5 = 1/2 (see FigEJ). 

The physics that we just described is very similar to 
the excitonic scenario for the resonance in the neutron 
scattering [j| and the pseudo-resonance in Bi g Raman 
intensity [t| in the cuprates. Like in cuprates, the res- 
onance in the two-leg ladders is the combination of the 
two effects: the presence of the gap in a single particle 
excitation spectrum, and the attractive residual interac- 
tion between quasiparticlcs. The attraction leads to a 
formation of a two-particle bound state below twice the 
gap, which shows up as a resonance in R(u>). The finite 
intrinsic width of the Raman resonance in a ladder is 
due to the fact that Raman intensity is a q = probe, 
and it includes a contribution from the states for which 
ImRo(uj) is finite below 2e m i n . 



In Sri4Cu2iC>4i, the sharp peak in the Raman pro- 
file has been detected around 3000cm" 1 0, 0. The 
measured maximum frequency 2e max is near 4000cm" 1 . 
Whether or not this peak is our pseudo-resonance de- 
pends on whether this material is truly described by a 
5 = 1/2 ladder, or Haldane effect is suppressed by 3D 
couplings, and the form of the excitation spectrum is sim- 
ilar to the quasiclassical expression. In the first case, nu- 
merical studies indicate [4[, |6j that the pseudo-resonance 
is too broad to account for the data. However, if the 
quasiclassical description is valid down to S — 1/2, the 
pseudo-resonance below 2e m ;„ = 2y/2 J1J2 is quite sharp 
(FigEH) and the profile of R(uj) is consistent with the 
data. Matching the peak position and the location of 
the upper edge for R(u>) by quasiclassical formulas yields 
Ji - 1600cm" 1 {200meV) and J 2 /Ji ~ 0.4 - 0.5. This 
value of Ji is comparable to J ~ lOOme!^ in the 2D 
cuprates [ijj, the value of Jij J\ roughly agrees with 
other estimates 0,0- 

In the rest of the paper we present the details of our 
derivation of the Raman intensity. Our point of depar- 
ture is the Hubbard Hamiltonian for two chains 

<i,j> i 

+ C/ E«T<l+<T<i)- ( 4 ) 

i 

In the quasiclassical case, we introduce antifcrromagnetic 
long-range order at Q = (ir,n) via < Y^k c k+Q,t c k,l >= 

a, and < J2k^k+Q t^M >— ~ a - Decoupling the Hub- 
bard U term using these relations, diagonalizing the re- 
sulting quadratic Hamiltonian, and introducing new va- 
lence (b, /) and conduction (a, e) electrons instead of 
c fc , c k+Ql d k ,d k+Ql we obtain 

H = J2 ^k(al a a kta -bl a b Ka ) + ^ k (el a e kia ~fl a f Ka ), 



where E fc = ^ (2t cos k + 1') 2 + A 2 E k = 

(2t cos k — t 1 ) 2 + A 2 , A = Ua, and prime indi- 
cates that the summation goes over magnetic Brillouin 
zone. The self-consistency equation on a yields a — 
f(t/U,t'/U), where /(0,0) = 1/2 0. 

The two-magnon Raman profile is observed in 
near-resonant Raman scattering regime, where the 
light couples to electrons predominantly via jA 
term, where A is the vector potential of light, 
and j is electron current with the components 



Jx = 



k.a 



2£sin(fc) 



a l,<r b k,a 



h.c 



h 



where 



k,a 

x and y are the directions along the chains and trans- 
verse to the chains, respectively. 

Other elements of two-magnon Raman scattering are 
the electron-magnon coupling, the magnon propagator, 
and the magnon-magnon interaction. They all are ob- 
tained in a straightforward manner from Eq. J3J| by tak- 
ing the large U limit, extending the Hubbard model to 
large S, and computing the spin susceptibilities in the 
RPA approximation, which becomes exact in the quasi- 
classical case [lfl, ll2|- The poles of the spin suscepti- 
bilities determine magnon dispersion, and the effective 
electron-magnon Hamiltonian is obtained by summing 
up RPA scries of particle-hole renormalizations of the 
Hubbard U. Once this interaction is known, one can 
straightforwardly obtain the effective vertex Mr{H, uj) for 
the interaction between light and two magnons with mo- 
menta k and —k and total frequency u> (see Ref. and 
Fig^J. Carrying out rather cumbersome calculations of 
the susceptibilities and vertices, we obtain at large U 



M R (k) = 



4Ji J 2 S'(cos(fc) - 1) 
24k) 



( G ? ; T G 



fx 



e iy e fy 



), (5) 



where e are unit vectors, we introduced J\ = -^ijj , Ji = 
4 g 2 ^- , and e(k) is given by Eq. OJ. 

The last element required for the computation of two- 
magnon Raman scattering is the magnon-magnon inter- 
action. To derive it, we note that at large U, the mag- 
netic properties of the Hubbard model are adequately 
described by the effective Heisenberg Hamiltonian H = 
Sr, M (JiSrSr+fi* + ^S,. S r+AIy ) . Applying the Holstein- 
Primakoff transformation to this Hamiltonian, and diag- 
onalizing the quadratic form, we reproduce the magnon 
dispersion and also obtain the four-magnon interaction 
vertex. In the quasiclassical approximation, only the 
term Hi nt = —V(k, l)e' k e_ k eie-i with two creation and 
two annihilation boson operators e is relevant [ill ]. For 
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FIG. 4: a) Diagrammatic representation of the Raman vertex 
for the interaction between light and magnons. Solid wavy 
lines represent light, solid and dashed straight lines represent 
conduction and valence electrons, and dashed wavy lines rep- 
resent magnons. There are other diagrams (not shown) whose 
role is to cancel parasitic contributions from these two dia- 
grams [Tl|. b) Diagrammatic derivation of Eq. Shaded 
rectangles represent V(k,l). 



V(k, I), we obtained 

V(k, = 1 i J 2 + Ji(cos(fc - I) + cos(fc + I))] 
+X 2 k X 2 ) + [J 2 + Ji(l + cos(fc - I))] fi k fiiX k \i 
— [3J 2 /2 + J 1 (cos(A;) +2cos(Z))] ^ k X k {tf + Xf) 



[3 J 2 /2 + Ji(cos(0 + 2 cos(fc))] w A, (fi 2 + X 2 ) (6) 



The coherence factors fi k and X k are given by 



1 2J 1 S + J 2 S 1 



1 2J 1 5'cos(fc) + J 2 S 2J X S + J 2 S 
k ~ V2\2JxScos(k) + J 2 S\y e(k) ' U 

With these results at hand, we can now compute the 
Raman intensity R(u>). Without magnon-magnon inter- 
action, the Raman intensity R(uj) = —ImRo(u>), where 



Ro(u) = du>'M 2 R {k)G k , u ,G- k ^ u ,, (8) 

and G kul = l/(uj—e(k)+iSsgnuj) is a magnon propagator. 
Substituting the result for M R (k) from J5J and evaluating 
the integral over uj' , we obtain 



R (lo) = (2J 1 J 2 S) 2 J2 



1 — cos(fc) 



uj - 2e(fc) + iS' 



(9) 

Like we said, the magnon dispersion e(fc), Eq. ijljl. has 
a maximum e max = 2J\S(1 + J 2 /2Ji) at k = ko = 
arccos(— J 2 /(2 Ji)), and a minimum e m i n = S\/8JiJ 2 
at k — 7r. Near the maximum and the minimum of 
e(k), the magnon density of states N(e) = dk/de(k) di- 
verges as a square-root of a distance to either e m i n or 



4 



(-max- Replacing J dk by J N(e)de, we reproduce Eq. 

Below 2e min , Rq(lj) is finite because of low-energy 
magnon states at small k, but is reduced due to the factor 
[(1 — cos/c)/e(fc)] 2 oc k 2 . We found from that as u> ap- 
proaches 2e m i n from below, ImRo(u)) tends to constant 
value. It then jumps to infinity at u> = 2e m i n + 0, and 
decays as 1 / \Jlo — 2e min at larger frequencies(see Fig. 01. 

When magnon-magnon interaction is included, the 
full R(u>) is given by Eq. J2J, where a — 
V(k,l)/(M R (k)M R (l))V 2 (see Fig. Hd). Near the two 
points where the density of states diverges, the inter- 
action V(k, I) from ©, (JTJ can be approximated by 
V(tt,tt) = ( Ji + J 2 /2)/2, and V(k ,k ) = J 2 (l + 
Ji/2J-\)/A. Using these forms, we obtain, e.g., near 

2£m.7.rj. • 



R(to) oc Im 



l + 2U(7r,7r) I(lo) 



1 



k~7T 



2e(k)' 
(10) 

and ui = u) + id. Comparing and (|l(Jf> . we see that 
there are two differences between Ro(uj) = —ImR$(uj) oc 
—Iml(uj) and R(uj). First, R(u>) vanishes at 2e m i n be- 
cause the divergence in ImI{oj) in the numerator of i|l(J|) 
is overcompensated by even stronger divergence in the 
denominator. By the same reason, R(u>) also vanishes at 
2e max (and at 2eq at finite S). Second, below 2e m i n , 
Iml(uj) is small and Rel(uj) is negative and behaves 
as — \/^/2e m in — oj. Then, at some w = u> res < 2e m i n , 
2V(ir,ir)ReI(u) res ) = 0, and R(lu) develops a pseudo- 
resonance. Near 2e max , Rel(ui) is positive, and the reso- 
nance does not develop. 

In Fig. we plot R(uj) obtained by solving Eq. (J2J) 
numerically using Eqs. @ and (Q) for magnon-magnon 
interaction, and Eq. (0) for the Raman vertex, and by 
formally extending the quasiclassical formulas to S = 
1/2. We see that R(uj) has a sharp peak slightly below 
2e m in- The position of the peak and its width somewhat 
depend on the ratio J%/J\, but for not very small J2/J1 
this dependence is rather moderate. For J2 = 0.5Ji, the 
peak is located near 2Ji, and its FWHM is 0.4Ji. In 
between 2e m i„ and 2e max , the intensity passes through 
a maximum, but R(u>) at the maximum is much smaller 
than at the peak below 2e m i n . Like we said, for S — 
1/2, numerical results indicate 0, [f| that the peak is 
more broad then in the quasiclassical analysis extended 
to S = 1/2. For larger S, however, the quasiclassical 
results should be more accurate. 

To conclude, in this paper we argued that the Raman 
intensity in spin S two-leg spin-ladder materials has a 
pseudo-resonance peak. The peak originates from the 
existence of a local minimum in the magnon excitation 
spectrum, and is located slightly below twice the magnon 
energy at the minimum. The physics leading to the peak 



is similar to the excitonic scenario for the neutron and 
Raman resonances in the superconducting state of the 
cuprates. At large S, the peak is quite narrow, its in- 
trinsic width scales as 1/S. For S — 1/2, though, the 
pseudo-resonance may be already rather broad. 
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